---
title: "PLS January Analysis"
output: html_notebook
---

The following replicates results for Ploger, Dunaway, Fournier and Soroka (2021), "The psychophysiological correlates of cognitive dissonance," forthcoming in Politics and the Life Sciences.

```{r setup}

library(stargazer, quietly = T)
load("PLS.ploger.workingdata.Rdata")

```


Videos are as follows:

11: Chimbolte (a town burns down)
12: May Day parades
13: famine in Niger
14: fighting in Sri Lanka
15: gorillas leave the zoo
16: a car that folds
17: kid makes movies
18: kid survives cancer


```{r Table 1}

model1 <- lm(v.gslmc.12 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model2 <- lm(v.gslmc.12 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model3 <- lm(v.rmssd.12 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model4 <- lm(v.rmssd.12 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
suppressWarnings(stargazer(model1, model2, model3, model4, no.space=T, type="text", omit.stat=c("f","ser")))

```


```{r Appendix Tables}

model1 <- lm(v.gslmc.11 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model2 <- lm(v.gslmc.13 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model3 <- lm(v.gslmc.14 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model4 <- lm(v.gslmc.15 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model5 <- lm(v.gslmc.16 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model6 <- lm(v.gslmc.17 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
model7 <- lm(v.gslmc.18 ~ female + age + educ + income + leftright, data=B[B$dropgsl.new==0,])
suppressWarnings(stargazer(model1, model2, model3, model4, model5, model6, model7, no.space=T, type="text", omit.stat=c("f","ser")))

model1 <- lm(v.gslmc.11 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model2 <- lm(v.gslmc.13 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model3 <- lm(v.gslmc.14 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model4 <- lm(v.gslmc.15 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model5 <- lm(v.gslmc.16 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model6 <- lm(v.gslmc.17 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
model7 <- lm(v.gslmc.18 ~ female + age + educ + income + wp.right2, data=B[B$dropgsl.new==0,])
suppressWarnings(stargazer(model1, model2, model3, model4, model5, model6, model7, no.space=T, type="text", omit.stat=c("f","ser")))

model1 <- lm(v.rmssd.11 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model2 <- lm(v.rmssd.13 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model3 <- lm(v.rmssd.14 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model4 <- lm(v.rmssd.15 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model5 <- lm(v.rmssd.16 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model6 <- lm(v.rmssd.17 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
model7 <- lm(v.rmssd.18 ~ female + age + educ + income + leftright, data=B[B$dropbmp==0,])
suppressWarnings(stargazer(model1, model2, model3, model4, model5, model6, model7, no.space=T, type="text", omit.stat=c("f","ser")))

model1 <- lm(v.rmssd.11 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
model2 <- lm(v.rmssd.13 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
model3 <- lm(v.rmssd.14 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
model4 <- lm(v.rmssd.15 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
model5 <- lm(v.rmssd.16 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
model6 <- lm(v.rmssd.17 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
model7 <- lm(v.rmssd.18 ~ female + age + educ + income + wp.right2, data=B[B$dropbmp==0,])
suppressWarnings(stargazer(model1, model2, model3, model4, model5, model6, model7, no.space=T, type="text", omit.stat=c("f","ser")))


```